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Using the effective action formalism, we investigate collective excitations of a dilute mixture of a 
Bose gas and a two component Fermi gas when both bosons and fermions have undergone superfluid 
transitions. We show that there is a repulsion between Bogoliubov and Anderson modes, which 
has important implications including disappearance of boson superfluidity. We derive an analytic 
expression for the long-wavelength dispersion relation of the mixture at zero temperature and give 
a condition for the instability. We also numerically calculate the damping rate due to boson- 
fermion interaction at finite temperatures and show that the two modes are stable at sufficiently 
low temperatures. 



I. INTRODUCTION 

Experimental realization of Bose-Einstein condensa- 
tion (BEC) of trapped atomic Bose gases and achieve- 
ment of trapped quantum degenerate Fermi gas (DFG) 
are two of the most important prospers in the field 
of ultracold atomic gas. Recently, '^''K atoms has been 
cooled to one-fifth of the Fermi temperature (0.2Tp) [§. 
One of the cardinal goals for the cooling of Fermi gas is 
the observation of BCS-type superfluidity . Except for 
^Li which has an enormously large negative scattering 
length [H , creation of effective attractive interaction for 
pairing is necessary, methods such as Feshbach resonance 
and phonon exchange in Fermi-Bose mixtures have 
been proposed 

Sympathetic cooling of atomic Bose-Fermi gas mix- 
tures has also risen to a thriving field. Since Bose gas 
can act as a coolant for fermions, severe decrease of 
rethermalization due to Pauli blocking can be avoided. 
In addition, this leads us to a fascinating system of di- 
lute Bose-Fermi mixture to explore. Success in simulta- 
neous trapping of Bose-Fermi isotopic mixture has been 
reported and simultaneous quantum degeneracy 

at temperature as low as 0.25Tf' has been observed in 
bosonic ^Li and fermionic ^Li mixture p^ . Density pro- 
files in harmonic traps |l3|Jl^ ] and collective excitations 
P, p"5p^ of Bose-Fermi mixtures has been investigated so 
far. 

Collective excitations of superfluid Fermi gas are pre- 
dicted to be detectable 0,|l^ and sound propagation in 
BEC has already been observed |]l9|. It would be in- 
teresting to investigate how these two sound modes act 
upon each other. In this paper, we study the collective 
excitations of dilute Bose-Fermi mixture when both bo- 
son and fermion have undergone superfluid transition. 
Spatial homogeneity is assumed throughout the calcula- 
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tion since overlapping of bosons and fermions can only be 
large in a gradual or box-like trap. Based on the effective 
action formalism, long-wavelength dispersion relation of 
the superfluid mixture at zero temperature is derived and 
damping rates due to coupling of the Bogoliubov and An- 
derson mode are also obtained by analytic continuation 
of the RPA-polarization bubbles. Instability of the su- 
perfluid modes in this mixture is also discussed. 



II. FORMULATION 
A. Model 

We shall start from a imaginary time path integral for 
which the grand canonical partition function reads 

Z = y I?0*2?02?V*2?'/'exp| - i(S'B[0*,(/)] 



where the total action consists of terms representing the 
Bose gas, 
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the Fermi gas, 
SF[iJ*,i'] = ^J dr J dx 



(3) 
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and boson-fermion interaction 

= V / dT f dx gMx,r)\^\^{x,T)\\ (4) 

where (/)(x. r) are complex fields and ?/'q,(x, r) arc Grass- 
man fields describing Bose component and the two hy- 
perfine Fermi component a = {t, i} respectively, gx = 
'^'^^^ is the coupling constant, for which ag, ap and | 
are the scattering length of boson-boson, fermion-fermion 
and boson-fermion interaction respectively. Note that 
only interaction between two di£Fc;rent hyperfine states 
is considered since Pauli principle prohibits s-wave scat- 
tering between two identical fermions. rrix = is 
the reduced mass and is the chemical potential of the 
corresponding component. 

In order to introduce an auxiliary field corresponding 
to BCS pairing field, we perform a Stratonovich- Hubbard 
transformation to the fermion-fermion interaction, 
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+A* (x, r) V'i {x, T)ij)^ {x, t) 



+V'*(x,T)V*(a;,r)A(a;,T) 



dx 



|A(x,r)|^ 
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Integrating over the fermionic field, we obtain an effective 
action, 
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+ SB[<p\4>]-fiTr[\n{-g-% (6) 



for which the inverse of Green's function in the third term 
can be expanded perturbatively. 
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xS{x — x')5{t — t'). 



Bogoliubov approximation is then applied to the Bose 
field and the pairing field is also separated into its average 
and fluctuation, i.e. 



'){x, t) = y/n^ + cj)'{x, t), 
Aa(x, r) = Ao -I- 5{x, r). 



(8) 
(9) 



Assuming quantum fluctuation to be small, the third 
term in (6) can be expanded, 

- fiTr[ln(-^-i)] = -fiTr[ln(-^o-')] + fif] -TrK^oS)^], (10) 



where the fluctuation part reads, 

1 fK + L + M 5{x,t) 



S(x,r) = - 



5{x,t) 



-K + L + M 



= --{Ka3 + {L + M)ao + 6{x,T)ai), (11) 



with, 



K = ij{x, r) 



mpv^ix, t) 
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/Ti^{<j)'{x,T)+<p'*{x,T)) 

(12) 
(13) 

M = gD [VHi(0'(x, r) + <P'*ix, r)) + 4>'{x, T)(f>'*{x, r))] , (14) 



+c^'{x,t)<P'*{x,t))], 
ih 

L = —[W ■ v{x,t) + v{x,t) ■ V], 



in which v{x,t) = 2mJ^^(^^'''^ ^^'^ = ^^{^^t) 

are the superfluid velocity and chemical potential for 
Cooper pairs respectively. The average of boson-fermion 
interaction, gA = ^^^^^ and the deviation of interaction 
strength from average g^ = il-^ has also been intro- 
duced in order to rewrite (11) in Pauli matrices. 

Since we are dealing with a homogeneous system, it 
is convenient to perform a Fourier transformation. The 
unperturbed Green's function reads. 



h / itujOn + et (fc) Ao 
D{k,uJn)\ Ao itkJn-^\.{k) 



(15) 



(7) with. 



B. Perturbation Expansion 

Since the trace of (7) is invariant under unitary trans- 
formation, we can perform a rotation in Nambu space, 
U{9) = e~5<^3S(a;,T)^ which corresponds to a gauge trans- 
formation. Separating the pairing field into phase and 
amplitude, i.e. A(x,r) = Aa{x,T)e^^p'^'^''^\ and choos- 
ing 6 = 6p leads to a new inverse of Green's function 
with a real energy gap. 



D{k,(Jn) = {ifvjJn - ^i{k)){ifiJ^n + (^^{k)) - |Ao|^, (16) 

where e„(fe) = - A«a + gans and a;„ = ^ is the 
even Matsubara frequency. 

Expanding (10) to the second order of fluctuation 
around its minimum leads to the dispersion relation of 
collective modes. The requirements of the linear terms 
of fluctuations to be zero {dS'^^^ = 0) are the well known 
Hugenholz-Pines relation 



HB = gsriB + gini + g^n^, 



(17) 
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with Oa{x,T;x,T) = Ua for the Bose field and BCS gap 
equation 



T{x, t; X, t) 



Ao 
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for the pairing field. 



C. RPA Polarization Bubbles 

Expanding (10) to the second order gives rise to various 
kinds of RPA Bubbles : 



go{k,uJn) 

ho{k,UJn) 



—— ^ T{k, UJn)^{p + fe, ), (19) 
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-^^^Ql{k,UJn)G]{p+k,U}^+U}^), (21) 



), (22) 



note that inverting Q and J- at (22) gives /cq (fc, w„). 
The remaining bubbles such as : 

^^Yl {P+^) Gl{k,^n)Gi{p + k,LJ,n+^n), (23) 

32(fe, w„) 

= ■^11 (P+|)'ei(fe,^n)ei(P + fc,t^m+t^„), (24) 

can be obtained as a combination of (19)-(22) by con- 
sidering small rotation U{x) = e^'^^^^^'' of phase, since 
gauge invariance gives rise to a Ward identity ||2^,|2l[ . 

The real part of the bubbles can be calculated in a 
purely analytical manner |^,|l| at absolute zero, the 
result of integration is independent of the difference — 
91 and " Mb 
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In order to obtain the imaginary part, we have to per- 
form an analytic continuation carefully by representing 
the Matsubara sum in the form of a contour integral . 
The final result is, 

p,rn 

= -j^^ j j delmGiR{p,e)lTsiFii{p + k,e + uj) 

(29) 



/ , he , h{e + Lu) 
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In the limit of long-wavelength, i.e. 



2kBT 



< 1 and 



^ 1, the imaginary parts of bubbles can be obtained. 
Assuming g-^ = gi and = for simplicity, we get 

N{0)ttAIlj 
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Im/io(fe, w, T) = Im/o(fc, UJ, T), 
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(31) 
(32) 



Im/o(fc, u, T) = .^Im(fc+(fc, uj, T) - k^ (fc, uj, T)), (33) 



where E = Ve^TAf and c = < 1. (33) can be 
verified from the Ward identity of o-q rotation j2^. The 
remaining integration is done numerically. 



D. Inverse Green's Function of Mixture 

The second order of fluctuation gives us the desired 
dispersion relation of the mixture. In analogy to the 
pairing field, we separate the Bose field to it's phase and 
amplitude as well, 

0'(fc,w) = 0A(fc,t^) + i0p(fc,w). (34) 

Introducing a vector notation, 

vl-t^ (^A(fc,c.), Mk,^), S{k,u;), 6{k,u:)), (35) 

the second order of fluctuation can be written in a matrix 
form, 

^^'^ = du^^G-'{k,uj,T)^, (36) 

k 

where is the inverse of Green's function of the mix- 
ture, 
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with matrix elements 
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The dispersion relation of the mixture can be found by 
requiring determinant of (37) equals zero. 



III. RESULTS AND DISCUSSIONS 

In the limit of long-wavelength linear dispersion, the 
real part of the dispersion relation in absolute zero can 
be derived analytically, 

1 f gsnB 
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' 1 / gsnB 
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g\N{Q)vlnB 
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(43) 



plus and minus sign in front of the square root corre- 
sponds to the two eigen modes of the mixture, Ander- 
son mode and Bogoliubov mode. We can see from the 
last term inside the square root that there is a repul- 
sion depending on average interaction between boson and 
fermion (^a)- If the Anderson velocity exceeds Bogoli- 
ubov velocity, this effect stiffens the Anderson mode and 
soften the Bogoliubov mode, and vice versa. A simple 
plot of gA dependence to this effect is plotted for a super- 
fluid mixture of fermionic ^Li (Fig.l) and bosonic *^Rb 
(Fig. 2) as an example. 

As the interaction becomes strong enough to satisfy 
the condition 



> 



47ras»Tijj 
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(44) 



the frequency of slower mode becomes a purely imagi- 
nary number. This corresponds to a instability, however, 
phase separation of boson and fermion occurs before this 
can be observed [||. 

Including the imaginary parts of bubbles, damping of 
the two modes due to boson-fermion coupling can be ob- 
tained by finding the numerical poles of (37), roots are 
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FIG. 1. Modification of Anderson mode in "^Li and "'Rb su- 
perfluid mixture varying the boson-fermion interaction. Den- 
sities of fermion and boson are taken to be 4 x lO^^cm"^ and 
lO^^cm"^ respectively, ao is tlie Bohr radius. 
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FIG. 2. Modification of Bogoliubov mode in "^Li and "'Kb 
superfluid mixture varying the boson-fermion interaction. 
Densities of fermion and boson are taken to be 4 x lO^^cm"'^ 
and lO^^cm"'^ respectively, ao is the Bohr radius. 



found to be in the form oi lo = (c— i^)k where the imag- 
inary part is k dependent. A plot of 7 in mixture of ^Li 
and *^Rb is shown in Fig. 3. Note that there is no damp- 
ing at absolute zero, the reason is that quasi-particle does 
not exist and therefore no real process of quasi-particle 
excitation can take place. The imaginary part is rather 
small even at relatively high temperature which implies 
that the two superfluid modes are stable. However, the 
effect of pair breaking and Landau damping |^ will be 
important near the superfluid transition temperature of 
fermion Tcf, which is out of the scope of our calculation. 
Bose-Fermi superfluid mixture near T^f will be discussed 
elsewhere. 
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FIG. 3. Imaginary part of Anderson velocity (7a) and Bo- 
goliubov velocity (7s) for — ai — 150ao in a ^Li - ^'^Rb 



superfluid mixture with density of 4 x lO^^cm 
respectively. 



and let'em" 



IV. CONCLUSIONS 

We have calculated the dispersion relation and damp- 
ing of a dilute mixture of boson and two component 
fermion. Repulsion between the Anderson mode and Bo- 
goliubov mode due to boson-fermion interaction is found. 
Instability of boson superfluidity is predicted in strong 
coupling regime, where phase separation is also predicted 

. Damping is found to be small in low temperature re- 
gion where pair breaking effects are negligible. 
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